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ABSTRACT 

Planetary rotation rates and obliquities provide information regarding the history of planet forma¬ 
tion, but have not yet been measured for evolved extrasolar planets. Here we investigate the theoret¬ 
ical and observational perspective of the Rossiter-McLauglin effect during secondary eclipse ( RMse ) 
ingress and egress for transiting exoplanets. Near secondary eclipse, when the planet passes behind 
the parent star, the star sequentially obscures light from the approaching and receding parts of the 
rotating planetary surface. The temporal block of light emerging from the approaching (blue-shifted) 
or receding (red- shifted) parts of the planet causes a temporal distortion in the planet’s spectral line 
profiles resulting in an anomaly in the planet’s radial velocity curve. We demonstrate that the shape 
and the ratio of the ingress-to-egress radial velocity amplitudes depends on the planetary rotational 
rate, axial tilt and impact factor (i.e. sky-projected planet spin-orbital alignment). In addition, line 
asymmetries originating from different layers in the atmosphere of the planet could provide informa¬ 
tion regarding zonal atmospheric winds and constraints on the hot spot shape for giant irradiated 
exoplanets. The effect is expected to be most-pronounced at near-infrared wavelengths, where the 
planet-to-star contrasts are large. We create synthetic near-infrared, high-dispersion spectroscopic 
data and demonstrate how the sky-projected spin axis orientation and equatorial velocity of the 
planet can be estimated. We conclude that the RMse effect could be a powerful method to measure 
exoplanet spins. 

Subject headings: planets and satellites: atmospheres, fundamental parameters — stars: planetary 
systems — techniques: spectroscopic - infrared: planetary systems 


1. INTRODUCTION 

Planetary rotation rates and axial tilts are critical 
parameters in determining seasonal climate v ariations 
(Williams & Kasting 1997 Cowan et al. 20121 and are 
relevant to th e planet formation and evolution history 
(Agnor 1999). By definition the planet rotational rate 
and axial tilt (e.g. obliquity) are respectively the time 
it takes for a complete revolution and the angle be¬ 
tween the planetary spin angular momentum and the 
planet’s orbital angular momentum. Planets with obliq¬ 
uity < 90° and > 90° are considered to have prograde 
and retrograde rotation, respectively. Planet spins are 
well-constrained for the Solar system planets, span ning 
a wide range of rot ational rates and axial tilts (Cox 
& Pilachowski 2000) and are considered to reflect the 
unique form ation a nd ev olutionary history of each planet 
(Laskar & Robutel|1993 ). Planets accumulate rotational 
angular momentum from the relative motions of accreted 
material. The stochastic nature of planetary accretion 
from planetesimals allows for a random component to 
the net spin angular momentum. Prograde angular mo¬ 
mentum (i.e. spin) could be accumulated by a planet on 
a circular orbit within a uniform surface density disk of 
small planetesimals. The planet clears a gap and thus 
accretes a larger fraction of material from the edges of 
its accretion zone. Retrograde spins are considered to 
originate from giant impacts during the early stages of 
planet formation. Therefore constraints on the planet 
spins are of high scientific interest, because deviations of 
the axes of rotation may have been caused by impacts of 
large bodies during their early history. Even the axes of 
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rotation of the gaseous planets may have been affected by 
impacts on their rocky cores before these planets accu¬ 
mulated their large atmospheres of hydrogen and helium. 
Such atmospheres are considered to be accreted hydro- 
dynamically, in flows quite different from those which 
govern the dyn amics of planetesimals an d lead to pro¬ 
grade rotation (de Pater & LissauerpOOl Faure| 2007). 

Exoplanet tidal theory predicts obliquity erosion, on 
time scales < 1 Gyr, for planets on < 10 day orbits 
around nor mal low-mass stars, i .e. preventing seasonal 
variations (|Heller et al. 2011a|b). Observational meth¬ 


ods have been proposed to probe the rotation rates and 
obliquities for exoplanets from oblateness measurements 


Seager 

(2002); Seager & Hui (2002); Barnes & Fortney 

(2003) 

Palle et al. (2008 

i; Kawahara & Eujn 

(2010); Fu- 

jii &; Kawahara |2(J12f, 

but require precisions in excess 

of ~ 0.1 ^mag. 

Currently a rotational 
one extrasolar planet - 

rate has been probed for only 
the young fast rotator /3Picb, 
le axial tilts have not yet been 
d high-dispersion spectroscopy 
near-infrared has recently be- 
iterising the atmospheres of hot 

(Snellen et al. 20141 whi 

measured. Ground-base 
(R > 20,000) in the 
come successful in charac 

Jupiters (Snellen et al. 

20101 Brogi et al. 2012; Birkby 

et al. 

2013 Brogi et al. 

2014 Schwarz et al. 

2015). At 

high spectral resolution molecular absorption bands are 


resolved into individual lines allowing their robust iden¬ 
tification by line matching with model templates. As the 
planet orbits its star, the radial component of the planet 
orbital velocity changes by tens of km s _1 , enabling a dis¬ 
crimination of the Doppler shifted planet spectrum from 
the steady telluric contamination. The planet signal is 
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then extracted by cross-correlating the data with model 
spectra obtained by mixing the expected spectroscopi¬ 
cally active trace gases in hot-Jupiter atmospheres and 
assuming a range of vertical temperature pressure pro¬ 
files. 


Kawahara (2012) considered the effect of planetary 
spin on the planetary radial velocity in dayside spec¬ 
tra of exoplanets, simulated the effect and concluded 
that planetary radial velocity could be a powerful means 
for constraining planet spins. In this paper we describe 
the potential of the Rossiter-McLauglin effect during sec¬ 
ondary eclipse ( RMse ) combined with near-infrared high- 
dispersion spectroscopy to provide constraints on exo¬ 
planet rotation and obliquity. 

The rest of the paper is organised as follows. Section [2] 
describes the rotational effect in transiting exoplanets. 
In Section [3] we derive the planet radial velocity anomaly 
due to the RMse e ffect. We choose the methodology of 


Ohta et al. (2005) who found an analytic solution for 
the RiVl effect at primary transits and show that their 
solution can equivalently be applied to the problem 
discussed here. Section [4] details the amplitude and 
shape of the RMse effect and discusses potential targets. 
Section [5] discusses the potential application of the 
method in light of the available and future instrumenta¬ 
tion. Finally Section [6] is devoted to our conclusions. 


2. RM EFFECT AT PLANET SECONDARY ECLIPSE 

In its nature, the RM effect arises becau se of stellar 
rotation (Rossitcr 1924; McLaughlin 1924). During a 
stellar occuitation or a planetary transit, portions of the 
rotating star surface are temporarily obscured, causing 
the removal of particular radial velocity components 
from the stellar broadened absorptio n lines leading t o 
a temporal radial velocity anomaly (Winn et al.J 2005). 
The same effect is expected to take place with~a rotating 
planet at secondary eclipse. When the planet passes 
behind the parent star the light from the approaching 
and receding parts of the rotating planetary surface 
sequentially enter/exit the geometric shadow of the star. 
The temporal block of light emerging from the approach¬ 
ing (blue-shifted) or receding (red-shifted) parts of the 
planet causes temporal distortion in the line profiles 
of the planet’s spectrum leading to an anomaly in the 
planet radial velocity curve. The shape and amplitude 
of this anomaly depend on the planet rotation rate, 
axial tilt and spatial orientation (i.e. planet spin-orbital 
alignment). Measurements of rotational rates and 
axial tilts are of high scientific interest as they provide 
constraints on exoplanet formation and evolution. 


3. FORMALISM OF THE EFFECT 

To describe quantitatively the radial velocity anomaly 
caused by the RMse effect we assume a two-body prob¬ 
lem with a central star and a planet of masses m s and 


m „ , respe ctively. We refer the reader to Fig. 1 in Ohta et 


al. (2005) for a schematic illustration of the top view of 


the planetary orbit and their equations (1) to (7). The 
orbital velocity of the planet as a function of time with 


respect to the star, up to O(e) as detailed in Murray & 


Table 1 

List of Notation 


Variables Definition Meaning 


m s 
m p 
P.. 
a... 
e ... 
uJ ... 
E.. 
n... 
M. 


i. 

f 

ri 



FTg-1' 
Fig-1" 
Fig.T 1 

Fig.T" 

Fig.T 1 

Fig.T' 
Fig.T" 
Fig.T a 
Fig. T a 
Fig. T a 


Orbital Parameters 

Star mass 

Planet mass 

Orbital period 

Semimajor axis 

Planet orbital eccentricity 

Argument of periastron 

Eccentric anomaly 

Mean motion 

Mean anomaly 

time of pericentre passage 

orbital inclination 

True anomaly 

planet to star distance 


I p . 

Up 

\p . 

R s 

R P 

V p . 


Internal Parameters of Star and Planet 


Fig. 1 a 
Fig. T a 
Fig.Tb 

If 


Si 

§3 


Planet spin-to-y-axis angle 
Planet angular velocity 
Sky-projected spin-orbit angle 
Stellar radius 
Planet radius 

Planet surface velocity, R P Q, P 


Mathematical notation 


Fig.|T|a, b Normal vector to the planet orbit 

§ 3 Position of the star 

§ 3 Position of the planet 

§3 Star to planet ratio R * / R p 

Eq See Fig. lc 

Eq F2 See Fig. T c 

EqT3 See Fig. Tc 

Eq FT See Fig. T c 


iu p . 
X s .. 

Xp.. 

7- 

Vs- 

xo-- 

20 " 

c- 


Note. — see Fig. 1 in Ohta et al. (2005)1 for all quantities 
marked with (*) and are detmed as expected. 


Dermott (1999) is 


frad p ~ --- na sin i [sin(M + uj) + e sin(2M + w)l, 

m s + m v 


with all quantities defined in Table 1. 


( 1 ) 


An eclipse or occuitation of a part of the rotating plan¬ 
etary surface causes a time-dependant asymmetry in the 
absorption/emission line profiles. These asymmetries re¬ 
sult in an apparent shift of the central spectral line po¬ 
sitions when the lines are unresolved. 

To describe the radial velocity anomaly caused by the 
planet rotation, similar tofOhta et al. (2005), we initially 
set the coordinate system at the star centre and its y-axis 
to coincide with the observer’s line of sight (Fig.[l]a) . The 
planet position is described with the coordinates (x p , z p ), 
corresponding to the orbit plane position and the planet 
impact parameter. 

For simplicity of the mathematical description of the 
problem, we choose a reference system (x ', z') centred on 
the planet and rotated such that the z'-axis is parallel 
to the rotation axis of the planet (i.e. parallel to f l p , see 
Fig. [I] b) and the rotation axis lies in the y'-zf plane. We 
define an angle X p between the sky-projected rotational 
angular velocity and the normal unit vector of the planet 
orb it, n n , s ee Fig.[lJ b. This differs from the definition of 
Ohta et al. (2005|), who assume A to be the angle between 
the sky-projected stellar rotation axis and the normal 





























































































A method to measure exoplanet spins 


3 


vector of the planetary orbit n p . 

In all calculations we ignore differential rotation of the 
planet surface as well as motions associated with atmo¬ 
spheric dynamics. A point on the surface of the rotating 
planet with coordinates (x', z') will move with a velocity 
v p given by 

v p = n p x' sin/p, ( 2 ) 

where is the angular velocity of the planet. The as¬ 
sociated radiation will exhibit a Doppler shift defined as 


Av Q p x r sin I p 
v c 


( 3 ) 


with respect to the observer along the y'- axis (i.e. the 
line of si ght). We re f er the reader to Sec. 3, Equations 


(14-19) in Ohta et al. (2005) for a derivation of the radial 


velocity profile for a star and adopt their expression ( 20 ) 
rewritten for the planet: 


A v p = — sin I p 


ff x'I(x',z')dx'dz' 
JJ I(x',z')dx'dz' 


( 4 ) 


Equation [4] relates the radial velocity change Av p and 
the line intensity I[x',z'). Fig. [2] illustrates the different 
cases of the RM effect during secondary eclipse. 

We evaluate the integrals assuming uniform model of 
the planet surface intensity I(x',z'). We ignore the role 
of planet limb-darkening as our goal is to estimate the 
first order rotational effect, and leave inclusion of the 
limb-darkening for future investigations. We also con¬ 
sider the star to be completely optically thick. 

At ingress and egress the position of the stellar disc 


satisfies the relation R s — R p < (xf + z' s 2 ) 1/2 < Rs + Rp. 
In order to simplify the computational task we rotate 
the coordinates in a time-dependent manner so that the 
stell ar centre is always located along the new x-axis, as 
in Ohta et al. (2005), see Fig.[l]b: 



1 ( < 
R P \Z x 's + z' s 2 \ Z ' S 



( 5 ) 


The position of the star is given by 



( 6 ) 


where 77 is defined as 


Vs = 



- 1 . 


( 7 ) 


The intensity on the uniform planet surface at (x, z) is 
given by 


I{x,z) = 


'Io, x 2 + z 2 < 1 and (x — 1 — y s ) 2 + z 2 > y 2 , 
^ 0 otherwise, 

where 7 = R s /R p , R p = 1. The moments of the intensity 
then are 


J(x', z ) dx' d z! = R 2 


nl 0 — / / I(x,z) dxdz 


, ( 9 ) 




Figure 1. (a): Schematic illustration of the planetary orbit plane, 
spin axis, and the observer’s line of sight; (b): Planet secondary 
eclipse ingress and egress phases and rotation axis; (c): A zoom of 
the planet and star configuration at ingress in the new coordinates, 
(see Table)!] for symbol definitions). 


and 


JJ x’l (x', z') dx'd z’ = — i JJ (x' s x—z' s z)I(x,z) dxdS. 

718 S ( 10 ) 

The range of the integrals is denoted with S and is de- 
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fined as the star-planet overlapping region (shaded area) 
in Fig.[l]c: 

/■l r\/l-X 2 rxo r\/l 2 -{x-Zs)' 2 

dxdz = I d x I d2+ / di / _ dz, 

IS Jx O j — y/l—x 2 Jx s —~f J — \/1 2 — (x— Xs) 2 

(H) 

The star and planet circles intersect at (xo,± 2 q), where 


x 0 = l- 


7 2 - Vs 


z 0 = \/l -xl = 


2(1 + Vs) 

Vh 2 -Vs)[{Vs + 2) 2 - 7^1 

2(1 + Vs) 


( 12 ) 

(13) 


These calcul ations are e quiva lent to the calculations pre¬ 
sented in Ohta et al. (2005) and to simplify the final 
result we introduce 


C = 1 + Vs - = 


2 Vs + 7 2 + V 2 


(14) 


2(1 dr Vs) 

Then equations [9] and 10 are analytically integrated as 


I(x, z)dx d2 = Iq 


and 


sin 1 zo—(l+Vs)zo+ r y 2 cos 1 C/7 

(15)" 


JJ (; xx' s — zz' s )I(x, z)dx dz = I 0 x' s (l +Vs) 


- ZoC 


+7 2 cos 1 C/'y 


(16) 
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<1 
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'■ 1 
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; t 




-2 0 2 -2 0 2 -2 0 2 
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Finally combining these two results we find the radial 
velocity anomaly of the planet (during ingress or egress) 
as a function of the star position {x' s ): 


A o , • T -zoC + 7 2 cos !(C/7) 

/\v n = il n x a sinl n -i-, 

7r-sin x z 0 + (1 +7y s )z 0 -7 2 cos -1 (C/ 7 ) 

(1 7) 

The fi nal result is equivalent to the result of Ohta et al. 

(2005) which is expected given the identical overlapping 
area between the planet and star circles. 


4. RMSE EFFECT AMPLITUDE AND SHAPE 

The shape and amplitude of the RMse effect are illus¬ 
trated in Fig. [2] for nine representative configurations of 
the planet spin-orbital alignments (i.e. impact parame¬ 
ter b = -fi-cosi and A p ), corresponding to prograde ro¬ 
tation. In case of retrograde rotation (i. e. X p > 90°) the 
radial velocity curves are inverted. The RMse computa¬ 
tion in Fig. [2] assumes a Jupiter-like planet, with v sin I p 
at the equator of 12.6kms _1 , on a 20 day orbit around 
a Sun-like star. We choose a period rather longer than 
the currently typical hot-Jupiter orbital periods because 
close-in exoplanets are expected to rapidly synchronise 
their rotation with the orbital period due to tides raised 
by the star on the planet. The time to spin-down the 
planet rotation is given by 



Figure 2. Illustration of planet radial velocity curve anomaly 
due to RMse effect for nine (representative) prograde spin-orbital 
alignments. The curves are plotted with a constant 10 kms —1 offset 
for clarity. Cases A to I correspond to the top planet-star config¬ 
urations. The curves flip and invert as X p increases to 360° and 
when b < 0. 


where Q pi R pi M p , w, ui s , M * and a p are the planet’s 
tidal dissipation factor, radius, mass, rotational angular 
velocity, and synchronous (or orbital) ang ular velocity 
(Goldreich & Soter 1966 Guillot et al.||1996 ). Assuming 
Q p ~ 1CF ( Correia fe Laskar|2010 ) and orbital periods of 
4 and 20 d, we find T syn ~ lCF and 10 9 yr, respectively 
implying that a ~ 20 d orbit required to have spin-down 
times comparable to the ages of the currently observed 
transiting hot Jupiters. 

Our choice for the value of v sin I p is driven by the 
correlation of the equatorial rotational velocities of solar 
system planets and their ma sses, suggesti ng that more 
massive planets rotate faster (Hug hes|2003 ). Currently a 
rotational rate has been constrained only tor the massive 
and young planet j3 Picb ( vsinl p ~ 25 kms-1), which 
is in line wi th the spin veloc ity-mass relation of the solar 
system ( Snellen et al.|2014 ). 

The maximum amplitude of the RMse effect is ex¬ 
pected for central eclipses (i. e. b = 0) and X p = 
0 or 180°. In the case considered in Fig. [2] the ampli¬ 
tude is ± ~ 6 km s' 1 . Assuming the v sin I p of /3 Pic b we 
find rotation effect with an amplitude of ~ 11 kms -1 . 
No effect is expected in cases where equal portions red- 
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Figure 3. Mock data (grey dots) illustrating the RMse effect, 
best-fit radial velocity curve (red lines) and radial velocity residu¬ 
als. Upper and lower pair panels illustrate the cases of a WASP-19 
system, assuming tidally synchronised prograde and a retrograde 
rotating planet, respectively. The blue symbols indicate the binned 
radial velocity curves by 2 min. 


and blue-shifted portions of the planet are eclipsed (i.e. 
b = 0 and A p = 90° or 180°). This is also the case where 
the planet rotation axis is normal to the plane of the sky. 

It should be pointed out that in all cases that include 
6 = 0, v sin I p is degenerate with X pi i.e. various com¬ 
binations of these two parameters can produce the same 
radial velocity amplitude. Importantly, it is the ampli¬ 
tude difference of the ingress and egress RMse, their signs 
and the shape of the radial velocity curve that constrain 
the planet obliquity and break this degeneracy. 

5. AN OBSERVATIONAL PERSPECTIVE FOR RMSE 


We produce mock data to illustrate the planetary 
RMse effect assuming an 8 and 40 m class telescopes, 
equipped with a near-infrared, high spectral resolution 
(i.e. R ~ 100,000) spectrometer. We assume that 
the planet orbit is constrained with enough precision 
to be subtracted prior to the search for the RMse ef¬ 
fect. We calibrate our simulation adopting a precision 
of ~ 5kms _1 as achieved by Birkby et al. (2013) at 
3.2 fim for water detection in the day side atmosphere of 
HD 189733b. We scale this precision by the square root 
of the number of collected spectra (i.e. 48), the ratio of 
the employed spectral coverage (assuming a hypothetical 
high-resolution spectrometer with a wavelength coverage 
of 400 mn) and a square root of the ratio of the target 
planet-to-star flux ratio and the planet-to-star flux ra¬ 
tio of HD 189733b, i.e. 1.3 x 10~ 3 from Birkby et al.| 
(2013). This factor plays a critical role as it accountsTor 
the strength of the planet signal, which is determined by 
the planet to star flux contrast. The contrast decreased 
from ~ 10 -3 to ~ 10 -6 when increasing the planet or¬ 
bital periods from ~ 0.5 day to 20 day. We also factor in 
the sampling rate with the telescope diameter and the 
brightness of the target host star. An account for the 
change of the planet flux during a secondary eclipse is 
also incorporated by factoring the radial velocity uncer- 
tainties by s quare root of the flux of a secondary eclipse 
Mandel & Agol (2002) model with unity out-of eclipse 
baseline and zero in-eclipse flux. 

Fig. [3] shows mock data generated assuming a 40 m- 
class telescope and a hypothetical planet host star of 
brightness K = 5.5 mag (i.e. similar to the currently 
brightest transiting hot-Jupiter host star HD 189733). 
Fig-! shows the results for a planet with the physi¬ 
cal properties of WASP-19 b, which is considered to be 
tidally locked. We assumed a prograde (A p = 35°) and a 
retrograde (A p = 215°) rotation with v sin J p = 8kms _1 . 
We perform a RMse fit to the mock data utilising a 
Levenberg-Marquardt le ast-sq ua res m inimization and es¬ 
timate X p and vsm(I p ), (jMarkwardt]|2009). We find the 
planet signal to be well-detected (~ 5-er) when combining 
nine secondary eclipses. 

We also explored the potential of an 40m-class tele¬ 
scope to detect the RMse effect in the currently bright¬ 
est transiting exoplanet host-star HD 189733 b. Assum¬ 
ing a synchronised planet rotation (i.e. RMse ampli¬ 
tude of ~ 1 kms -1 ) and A p > 0° we find that ~ 50 sec¬ 
ondary eclipses need to be co-added to detect the spin of 
HD 189733 b. 

Finally, we investigated the potential of an 8m-class 
telescope to detect the RMse effect in K=5.5 host star 
with a WASP-19b-like planet and synchronised rotation. 
We find a number of ~ 50 secondary eclipses will have 
to be added in order to detect the RMse effect at ~ 3-cr 
confidence. 

6. DISCUSSION AND CONCLUSIONS 

We have discussed the Doppler signature of a spinning 
exoplanet during secondary eclipse and derived a rela¬ 
tion for the radial velocity as a function of the planet’s 
rotational rate and axial tilt (Eq.[l7|). The RMse effect 
complements the tool-box offered Toy transiting exoplan¬ 
ets providing a proxy to exoplanet spins. We showed that 
the formalis m of the RM se is eq uivalent to the formalism 
developed by Ohta et al. (2005) to describe the RM effect 
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caused by a planet transiting its star jus t from a different 
perspective. We note that Hirano et al. (2011) derived an 
improved accuracy radial velocity anomaly curves for the 
R M effect duri ng pri mary transits. However, the solution 
offOhta et al.| (2005]) provides a conservative estimate of 
the radial velocity amplitude and is precise enough to 
illustrate the RMse effect. 

We acknowledge that our assumption of a rigid ro¬ 
tating non-limb darkened planetary surface (i.e. ignor¬ 
ing differential rotation and atmospheric dynamics) is 
a crude approximation, but it illustrates the expected 
RMse radial velocity anomaly in the simplest case. The 
radial velocity originating from the RMse effect is ex¬ 
pected to be degenerate with atmospheric dynamics. Ex¬ 
oplanets are expected to have atmospheric circulation 
with in creasing dynamics at shorter orbital distances 
(Kataria et al. 2015). Therefore it is expected that the 
radial velocity curves of most close-in exoplanets may 
significantly differ from that assumed here, with super- 
rotating winds expected to enhance the radial velocity 
signals by factor of a few. The RMse effect presents 
the opportunity to probe atmospheric winds with alti¬ 
tude, when wide coverage high dispersion spectroscopy 
becomes capable of providing per-point precisions in the 
~ms' domain. In such cases it could be possible to 
measure longitudinally-integrated radial velocities and 
map the surface of a planet during the secondary eclipse 
ingress/egress phases, providing constraints on the hot¬ 
spot position. We postpone the derivation of the RMse 
effect in these cases for future theoretical global circula¬ 
tion model investigations and identify the top three plan¬ 
ets expected to exhibit the strongest RMse effect assum¬ 
ing synchronised rotation. Applying Eq. 17 (assuming 
A p = 0°) we find values of Au p ~ 3.7, 3.4, and 3kms _1 
for targets WASP-103 b, WASP-12 b and WASP-19 b, re¬ 
spectively. 


Kawahara (2012) demonstrated the rotational signa- 
ture of an exoplanet in the planet radial velocity. The 
RMse effect discussed in this study is of similar nature. 
However the geometric shadow of the star rather than 
the changing planet phase is the factor removing ve¬ 
locity components from the visible planet surface. Im¬ 
portantly, because the RMse effect occurs only during 
ingress and egress, its detection requir es significan tly less 
observing time compared to the case in Kawahara (2012). 
RMse therefore could be much-well suited for planets 
with longer orbital periods. 

We also demonstrated that to be put into practice the 
RMse effect requires exoplanet host stars brighter than 
the currently known exoplanet hosts, i.e. brighter than 
K ~ 6 mag and large aperture telescopes (ie. ~40m). 
The strongest RMse effect for evolved stars is expected 
for non-synchronized rotation, i.e. planets on longer than 
20 day orbits, with a typical z^sin I p 12.5 km s -1 or more. 
However, such planets exhibit lower temperatures (i.e. < 
700K), compared to the typical hot Jupiters giving a 
small planet-to-star flux contrasts (i.e. ~ 10 -6 or lower 
at ~ 3/xm). The detection of the RMse effect in such 
cases would be an extremely challenging if not impossible 
task, even for the upcoming 40 m class telescopes. 

An interesting case could be young planets which still 
have not been synchronised. Observations using high- 
dispersion, near-infrared spectroscopy have currently 


constrained the rotational rate of the only one such case 
(0 Pic b) from rotat ionally broadened absorption lines 
(Snellen et al. 2014). An important opportunity could 
exist if the orbit orientation of this planet allows transits 
and secondary eclipses, because the large brightness of 
the host star (i. e. K ~ 3.5 mag) and the fast spin of the 
planet (i. e. v sin I p ~ 27kms _1 ) would produce a strong 
and detectable RMse signal. Another interesting oppor¬ 
tunity could be offered by transiting brown dwarfs with 
orbital periods larger than 10-20 days. Such systems as¬ 
sume high flux contrasts between the star and the brown 
dwarf in the near-infrared which in case of a bright host 
star (i. e. K ~ 5 mag) could provide an opportunity to 
probe the spins and latitudinal radial velocity maps of 
these objects. 

In the near future, the Next Generation of Transit Sur¬ 
veys (NGTS), the Transiting Survey Satellite (TESS) 
and the Planetary Transits and Oscillations of stars 
(PLATO) projects are expected to significantly expand 
the sample of known planets hosted by bright stars (i.e. 
brighter than K ~ 9.5 mag) and hence to provide more 
targets suitable for detection of the RMse effect. 

The hypothetical radial velocity data in Section[5] as¬ 
sumes retrieval of the planet signal via cr oss-correlatio n 
with a model planet spe c trum similar to ISnellen et al7 


d2010|);|B r ogi et alT] (2012 ; Birkby et al.|Q2013P;|de Kok et 
aT ( |2U13| ); |Schwarz et al. (2015). Although currently the 


cross-correlation technique has not been demonstrated 
to provide per spectrum r adial velocity measurements, 


in their Fig. 1, Snellen et al. (2010) presented results 
for the planet geocentric radial velocity as a function 
of the planet orbital phase at a high significance (~ 3- 
er). It is expected that large-aperture near-future tele¬ 
scopes (i.e. 30-40 m) could provide even higher signifi¬ 
cance signals and planet radial velocity measurements 
from multiple spectra. In addition, cloud-free exoplanet 
atmospheres are expected to result in stronger cross¬ 
correlation signals compared to hazy or cloudy atmo¬ 
spheres for the same stellar and planetary physical prop¬ 
erties. Cloud-free atmospheres are expected to have well- 
pronounced lines, i.e. non-muted absorption features 
which are expected for hazy atmospheres. 

Stellar activity and pulsations are known to produce 
radial velocity variations that could in some cases mimic 
those induced by the orbital motion of exoplanets. This 
may lead to misinterpretations of radial velocity varia¬ 
tions, especially when those variations h ave peri o ds less 
or eq ual to the star rotational periods (Lagrange et al. 
2010). Therefore, it is important to consider whether 


stellar activity, for instance the distortion to the stellar 
spectral lines induced by starspots on a rotating stellar 
surface, could induce residual signal misinterpreted as ex¬ 
oplanet signal. This is especially pertinent for co-added 
observations obtained over several nights. 

A typical RMse observation would last for ~ 0.8 h on 
each of the ingress and egress phases of a WASP-19b- 
like planet (see Fig. [3]). This time interval is significantly 
shorter than the rotation periods of t he t yp ical p lanet 
ho st star s ranging from 10 to 40 days (Paz-Chinchon et 


al. 2015). Furthermore, starspots are expected to pfcT 
duce sine-wave radial velocity variations (i.e. different 
from the RMse shape) with amplitudes of a few hundred 
ms -1 in the optical and are expected to reduce by fac- 
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tor of at least a few at near infrared wavelengths, where 
the contrast between the photosphere and cool star spots 
(e.g. AT ~ 550K) is significantly reduced. 

Thus, if present, residual stellar signal caused by 
starspots would be a factor of a few tens smaller than the 
the expected 1-2 km s -1 amplitude of the RMse effect for 
the closest synchronized planets and even more for non- 
synchronized planets on longer than 20 day period orbits. 
Only in the case of most evolved planets (i.e. typical hot 
Jupiters), if synchronized, with orbital/rotational peri¬ 
ods from 4-5 days to < 20 day could a residual introduce 
larger radial velocity scatter when multiple observations 
are combined. However, such planets are expected to 
have RMse amplitudes much smaller than lkms -1 and 
would be difficult targets by definition. 

In the near infrared the planet spectrum is dominated 
by thermal radiation and is not expected to contain infor¬ 
mation from distortions caused by starspots. In case that 
the planet spectrum also contains a non-negligible com¬ 
ponent from a reflected star spectrum then the planet 
spectral lines will have distortions caused by the star 
spots. This would make the planet spectrum different 
and more difficult to cross correlate with model spectra. 
However, such an effect is not expected to mimic the 
RMse. 

In conclusion, we have demonstrated the RMse effect 
and its potential to constrain planet rotational rates and 
sky projected spin-orbital alignments. We derived the 
radial velocity curve caused by the RM anomaly and es¬ 
timated the amplitude of the effect for a hot Jupiter. 
Finally we discussed the prospects for detecting the ef¬ 
fect and constraining planet spins and axial tilts from 
the current and upcoming instrumentation. 
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